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O 

fvi , In a tunneling experiment across a quantum dot it is possible to change the coupling between the dot 

and the contacts at will, by properly tuning the transparency of the barriers and the temperature. Gate 
voltages allow for changes of the relative position of the dot addition energies and the Fermi level of the 
leads. Here we discuss the two limiting cases: weak and strong coupling in the tunneling Hamiltonian. In 
the latter case Kondo resonant conductance can emerge at low temperature in a Coulomb blockade valley. 
We give a pedagogical approach to the single channel Kondo physics at equilibrium and review the Nozieres 
scattering picture of the correlated fixed point. We emphasize the effect of an applied magnetic field and 
D ' show how an orbital Kondo effect can take place in vertical quantum dots tuned both to an even and 

5_i , an odd number of electrons at a level crossing. We extend the approach to the two channel overscreened 

Kondo case and discuss recent proposals for detecting the non Fermi Liquid fixed point which could be 
reached at strong coupling. 
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Recently systems have been fabricated which can sustain quantum coherence because of the smallness of their 
size, provided that the temperature is low enough. These mesoscopic systems are nanostructured devices in which 
quantum coherence sets in at very low temperatures and modifies the properties of the device as a whole. This happens 
notwithstanding the fact that the system is connected to an external biasing circuit. They have global, geometry 
dependent properties and striking quantization phenomena arise which are largely independent of the specific sample 

VO ' involved: charge quantization (in unit of the electron charge e), periodicity in the magnetic flux quantum 4> = hc/e, 

< ^-P conductance quantization (in units of Gk = 2e 2 /Trh = (6.5i"ff2) _1 ). 

Among these macroscopic quantum phenomena there is the unitary limit of the Kondo conductance in tunneling 
across a quantum dot (QD) at Coulomb Blockade (CB) [1] [2] that was first measured in 1998 [3]. The Kondo 
phenomenon is well known since the sixties and explains an anomaly in the temperature dependence of the resistivity 
of diluted magnetic alloys [4] [5] . 

This review is devoted to some features of equilibrium Kondo conductance across a quantum dot in the CB regime 

interacting with two contacts. It is remarkable here that the dot acts as a single impurity probed by the metal 

contacts, so that the properties of the Kondo state are not mediated over many impurities per cubic centimeter as 

it happens in diluted magnetic alloys. A strongly coupled state sets in between the dot and the contacts and phase 

JL" ' coherence is established among the conduction electrons and the whole structure. 

. £h ! The temperature scale for the interactions between a magnetic impurity and the delocalized conduction electrons 

of the host metal is the so called Kondo temperature Tk- It is defined as the temperature at which the perturba- 
tive analysis breaks down. Therefore, different approaches are required to investigate the thermodynamics and the 
transport properties of a quantum dot in the Kondo regime in the whole range of temperatures from the perturbative 
region down to the unitarity limit. 

Recently, the Kondo model and the Anderson impurity model in its Kondo limit have been deeply investigated 
by numerical renormalization group (NRG) calculations [6], the Bethe ansatz method [7] and conformal field theory 
(CFT) techniques [8]. Further developments in the NRG methods have been applied successfully in the crossover 
region T ss Tk [9] . The zero field spectral and transport properties of the Anderson model in the Kondo regime [10] 
as well as the field and temperature dependence of the Kondo resonance and the equilibrium magnetoconductance 
of the dot [11] have been investigated. The tunneling conductance as a function of the gate voltage has also been 
calculated with NRG, in wide temperature range for a single quantum dot with Coulomb interactions, assuming that 
two orbitals were active for the tunneling process [12]. 
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We leave out the case when the electron distribution is not in local equilibrium about the Kondo impurity and the 
linear response theory is no longer sufficient. A number of techniques have been applied to describe nonequilibrium 
properties such as the nonlinear conductance, the nonequilibrium stationary state or the full time development of 
an initially out-of-equilibrium system: variational calculations [13], perturbation theory [14], equation of motion [15], 
perturbative functional integral methods [16], non-crossing approximation (NCA) [17] [18], perturbative renormal- 
ization group (RG) method in real time [19]. In particular, the last technique is well suited to describe quantum 
fluctuations which are induced by strong coupling between a small quantum system and an environment. It succeeds 
in reproducing the anomalous line shapes of the conductance observed in several recent experiments [3] , due to the 
renormalization of the resistance and the local energy excitations [20] . 

We also leave out situations in which the levels localized at the dot are close in energy to the Fermi energy of the 
contacts (mixed valence models). 

The plan of the paper is the following. We start from the tunneling hamiltonian formalism when the coupling 
between the dot and the contacts is weak (Section II). A scattering approach in one-dimension is particularly 
suitable when studying the linear conductance in the device. The prototype model to account for strong electron- 
electron repulsion on the dot is the Anderson Hamiltonian with onsite repulsion (Section III). In the limit of strong 
correlation between the dot and the leads the latter model maps onto the so called "Kondo Hamiltonian model" 
(Section IV). A poor man's scaling approach up to second order leads us to the definition of the Kondo temperature 
Tk- Next, we introduce the physics of the single channel (Section VI) and the two channel Kondo problem (Section 
VII), both in the perturbative region and in the unitarity limit. In such a context we briefly sketch the Anderson, 
Yuval and Hamann Coulomb Gas approach [21] to the isotropic and anisotropic one channel Kondo system which 
gives a straightforward, though qualitative, insight of the strongly correlated state. 

In the conventional setting the dynamical variable which is coupled to the electrons propagating from the leads is 
the total dot spin. There are cases in which the spin is locked to orbital degrees of freedom or even absent ("orbital 
Kondo") [22]; such cases are better found in a vertical geometry in presence of a magnetic field orthogonal to the dot 
[23]. A rather strong vertical magnetic field can induce level crossings in the dot due to orbital effects and produce 
accidental degeneracies which can give rise to exotic Kondo couplings (see Section VLB — C). Some attention is 
drawn to the cylindrical geometry and to symmetry selection rules in the cotunneling process, also in connection 
with proposals for achieving the two channel Kondo non Fermi liquid fixed point in a vertical dot device [24] [25] [26] 
(Section VII.C). 

II. TUNNELING HAMILTONIAN AT WEAK COUPLING 

Quantum dots (QD) are fabricated in semiconductor hetcrostructures, by applying metallic gates to confine few 
electrons [27]. Because the confining area is quite small (the radius is ~ 100 -j- 100CL4), the charging energy is much 
larger than the energy associated with thermal fluctuations, provided temperature is below 1°K. Therefore the dot 
is only weakly linked to metal contacts and one can bias the system in such a way that the electron number N in 
the dot can be changed at will. Because of the confining potential, dots display a level structure organized in shells, 
exactly the same as atoms do. Analyzing this level structure is of primary interest because these few interacting 
electrons (N < 20) confined in a disk-like box (see Fig. 1), at special values of the parameters, arc ruled by strong 
many-body Coulomb correlations. Hund's rules can be seen at work: the total spin S of the electrons confined in 
the dot is maximized as long as no magnetic field is applied. On the other hand, a magnetic field in the direction 
orthogonal to the disk produces strong orbital effects which favour larger values of the total angular momentum M 
as well as total spin S. 

Quantum dots are remarkable because of Coulomb oscillations. At very low temperature linear conductance (at 
vanishing bias V sc i) is zero, except for peaks at discrete values of the gate voltage V g , when it is energetically favourable 
to add one extra electron to the dot. Therefore V g controls the number of particles on the dot in the CB regime. 

Wc tunc the chemical potential // of the left (L) and right (R) bulk contacts within the CB valley of the conductance 
at N particles, that is /j,n < /j, < HN+i,a = N+1 E a — N Eo; here the energies N E a are the total energies for the 
dot with N ± 1 particles characterized by the quantum numbers a and n Eq is the ground state energy (GS) at N 
particles. If N is odd, then the GS is at least doubly degenerate because the odd particle state can have spin a =f, {. 
In this case the dot can be treated as a single Anderson impurity. 

The peculiarity of these artificial atoms is in that the level structure can be investigated by measuring the tunneling 
current. Peaks in the linear conductance versus gate voltage V g separate regions at Coulomb blockade (CB) in which 
N is fixed and differs by just one electron. The peaks at zero source-drain voltage occur whenever V g compensates 
the required chemical potential for electron addition or subtraction. Hence, measuring a tunnel current across the 



device vs V g provides the spectroscopy of the dot levels. Adding a magnetic field the spin state of the dot can be 
changed, which in turn changes selection rules for electron tunneling [28]. 

The height and the width of the conductance peaks at resonance depends on temperature T. Let the dimensionlcss 
conductance be g (in units of 2%- (factor of 2 is for the spin)). In the simplest case the maximum of the conductance 
at the resonant peak is g m ax oc r/(fcsT) and the halfwidth of the peak T — ttv (0) \V\ is ex T; here V is the 
tunneling strength (see below) and v(0) is the density of states at the Fermi energy. In the CB region tunneling via 
virtual states of the dot with N ± 1 is only fourth order perturbation in V (they are named cotunncling processes), 

v '' — — 1 /h, and vanishes exponentially with reducing temperature. 

Nevertheless, unexpectedly, transport measurement in a dot gives rise to some new physics which can be related to 
the Kondo effect of spin impurities in non magnetic metal alloys (see Fig. 2) [3]. 

In the following we introduce the model for a quantum dot interacting with the contacts and define an equivalent 
one-dimensional Hamiltonian. 

A. Current within the tunneling Hamiltonian 

In this Subsection we discuss the mutual influence between the contacts and the dot, starting from the weak link 
limit. Conduction electrons in the leads constitute a Fermi sea (FS) of non interacting fermions with plane wave single 
particle wavefunctions on side L and on side R. The dot is described by an Hamiltonian Hjj and coupling between 
the dot and the leads is accounted for through a tunneling term, so that the Hamiltonian describing the whole system 

is: 

H = H D +Y, e kA* a k* + J2 e kA a b k « + ^[V ka c\ aa da + V^dUkaa}- (1) 

ka ka kaa 

The right and left Fermi sea (FS) of the two contacts R, L are acted on by operators c kRa = a ka and c k L a = b klJ ( 
here the index a stands for R, L) . The canonical transformation of the dot problem made by Glazman and Raikh [1] 
is just the construction of two species of fermions whose wavefunction is even or odd with respect to the dot center, 
in the case that the two barriers are unequal V a = V r ,Vl. It changes the picture from operators a ka and b ka to 
operators a k(T and fi ka given by 

a klJ = ua klJ + vb k(7 , flko- = ub ka - va klJ 
V R V L 



u=-f, v=^, V=^\V R \2 + \V L \i, T RyL = ir\V R , L \ 2 v(0) (2) 

where i/(0) is the density of states at the Fermi energy. 

A general formula for the the conductance of interacting systems [30] [15] [31] has been written, resorting to the 
nonequilibrium Keldysh formalism [32] ; in such a framework the current through the interacting region is written in 
terms of the distribution functions in the leads and local properties of the intermediate region, such as the occupation 
and the density of states at the dot site: 

J=^Jduj(Tr {2 [f L (u)T L - f R {w)T R ] (G r - G a )} + Tr{2(T L - T R ) G< }) , (3) 

where G r , G a , G < are the usual retarded, advanced, Keldysh Green functions for the dot in interaction with the 
leads and the / are the Fermi functions. The Green's function G r will be denoted as Gdd in the following. Both G r 
and T R l are matrices in case that many channels are present. 

Formula (3) can be cast in a more simple form in the case that the couplings to the leads differ only by a constant 
factor [15]: 

^ fdoj [f L (u>) - Jr{uj)] ^mTr {fG»} , (4) 
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where f = Y R Y L /T. Here r = T R + Y L . 

Now we derive explicitely Gdd within the lowest order perturbation in the tunneling. 

We start from the Hamiltonian in eq. (1), where Hjj is given by a single impurity energy level td- All quantities 
will be scalar quantities for simplicity. We define the electron retarded Green's function for the unperturbed leads: 
Gq" 1 = J2 k (uJ — £fe + i0 + ). Projecting the equations for the total Green's function: 



(iw n - H)G(iu> n ) = 1 

G\iLu n ){-iw n -H) = l (5) 

onto states in which one single extra particle is occupying the delocalized state \k > or the impurity state \d > we 
have: 

(iu> n - e k )G ktk > a (iw n ) = Skk 1 + VkGd,k'cr(iv n ) 

Gd,k'<r(iWn)(-iU n ~ Cfe') = G d ,d<r{iuJn)V k * (6) 

where G d ,ka(iu) n ) =< d\G\k >, G ktka (iu> n ) =< k\G\k >, G d ,da{iu n ) =< d\G\d > and V k = (k\V\d). 

This gives for each scattering channel (here in order to simplify the notation we suppress the channel label a): 



= G ka 5 kk , + G ka VkGd,daV k *G k , a ,. (7) 

The density of states of the scattering electrons is defined as: 

K^) = --^E G UH- (8) 

77- t J > 



IT 

k 



Using similar equations for the dot, we write: 

Gd,da{iUn) = -. ^ i T/ i 2 //■ r- (9) 

iLu n - e d - 2^ fe \Vk\ 2 /(iuj n - e k ) 

Now, let us consider the L and R contacts as two equal Fermi gases of noninteracting particles and equal chemical 
potential /i = (within linear conductance regime a source drain voltage V s d is not applied). In such a case it is 
enough to discuss a single effective contact and the corresponding wavefunctions are plane waves of wavevector k 
in the x-direction with a label a which includes all other quantum numbers. Their energy dispersion e ka can be 
linearized about fj, with e q w hvpq, where q is the momentum measured with respect to the Fermi momentum and 
vf is the Fermi velocity. Using a constant density of states i/(0) (at the Fermi energy per spin L /2irhvF where L 
is the size of the system) and a bandwidth of size 2D and neglecting the k dependence of V k , the sum in eq. (9) is 
readily done: 

} -. = \V\ 1/(0) / de- = -ln- = S d («w„). (10) 

^ iuj n - e k J_ D iu n - e n iu n - D 

The retarded Green's function for the dot is obtained from Gd,da{iu n ) in the limit to real frequencies: iuj n — » lu + «0 + 
according to: 



G d,da(* w ") = {* w « - e<i - Sd(iw„)} -» cj - e d In 



v^?D + w 



vpD — to 



iY = w-e d + iT, (11) 



where td is the renormalized dot energy The coupling of the dot to the contacts shifts the location of the pole 
corresponding to the energy of the localized level and gives a finite width Y to the resonance. 

B. 1-d Scattering formalism 

At zero temperature scattering is only elastic. In the following we develop a 1-d scattering approach which 
is mostly useful in vertical structures [33], next we show that also the tunneling conduction can be cast into the 
scattering formalism. 

If the evolution operator U(t,t') is known, the transmission amplitude can be written in a scattering approach as 
9 R ^ L =< b 1f U(+oo,—oo)a >= v*u < aJ(oo) a(— oo) > —u*v < f3' (oo) /?(— oo) >= v*uS° + u*vS x , where we have 
defined the two expectation values as the scattering matrix elements for the two uncoupled even and odd channels. 
Hence, the conductance takes the form of the Landaucr formula: 



T 



R^L _ \ a R->L\ 



A\v*uY 



^ 



4r fl r £ 



&* 



(12) 



The potential after the transformation of eq. (2) has become even, hence unitarity is satisfied by each channel 



(Tr+Fl) 



-. In particular, if the potential 



individually: \S l \ 2 = 1, / = 0,1. The unitary limit of the conductance is g u 

is S(x), odd parity is totally transmitted and such that S 1 = — 1. 

Eq. (12) is valid for a general system with an interacting intermediate region. One can match this approach to a 
1-d scattering approach for non interacting electrons as it is done here below. 

Let us place the impurity (QD) at x = and consider the scattering amplitude Jl,r of a one electron wavefunction 

ip(x) 



4>>(x) oc e 



tkx 



■ fne 



ikx 



^ < {x)^e lkx + .f L e 



— ika. 



x » 
x « 0. 



(13) 



Here the transmission coefficient is T = |1 + /r| while the reflection coefficient is R — |/l| 2 and they satisfy the 
conservation of flux: T + R = 1. 

If the dot structure is even with respect to the origin along the vertical axis, the even parity I = and the odd parity 
I = 1 channels are independent. It is then useful to define even and odd parities f l : f° — |(/l + /_r), f 1 — \{Jr — /l) 
and the elastic t— matrix t l — ik f l /tt. Here the energy of the incoming particle is hvpk (in units Tivf = 1, it follows 
that [t] = energy) and k = 2n/L (where L is the linear size of the system). The t— matrix is related to the S— 
matrix according to: 



5' = e 2lS 



S l - 1 = - 



2m . 



■f 



f = 



— sva5 L e w 



(14) 



where 8 are the phase shifts for the two parities. In this context unitarity of the S— matrix, | J^ IS*'! 2 = 1, is the 
same as flux conservation R + T = 1 . 

Conductance is given by the Landauer formula: 



conductance = g = T = \l + fn\ = 



1 + i^^sinS 



l e iS 



(15) 



In the case of resonant tunneling we have T ~ 1 and R = \Jl\ 2 ~ 0, so that eq. (15) becomes g — > 1, that is the 
unitary limit. 

The condition R = \n(t° - ^/(iko)] 2 = sin 2 ((5° - 5 1 ) ~ yields (5° - S 1 = Smod tt, while in eq. (15) wc have T ~ 1 
for (5° = S 1 = 6 -» tt/2. 

If the potential is even (r^ = T^), unitarity is satisfied by each channel individually: \S l \ 2 — 1, / = 0,1; in 
particular, if the potential is S(x), odd parity is totally transmitted and such that S 1 = — 1, and the conductance 
becomes: 
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i<s"-i: 



— f 



sin 2 5°. 



(16) 



Again, if 5° = 7r/2 conductance reaches the unitarity limit. 

In the following we describe the basic approximations and the weak coupling limit of very opaque barriers between 
dot and contacts. 

In this case tunneling of lead electrons across the dot is a perturbative process as derived in Section IIA. Within 
such a perturbative regime neither the dot nor the contacts are much affected by the interaction with the other. The 
relevant effect on the QD is the shifting of its levels and a level broadening, which is a second-order effect in the 
tunneling strength V. This shows itself as a small imaginary part added to the energy levels, i.e. a finite lifetime. 

It is easy to show that the linear conductance given above can be rephrased in terms of the imaginary part of the 
Green's function as in eq. (4) at the perturbative level. To this end, let us first state some formalities regarding the 
selfenergy S and the t matrix: 

(E-H - E)G = 1; G (£ - H ) = 1 -» GoG" 1 + G S = 1 -> G = G + G SG. 



In the Born approximation X and t coincide because by definition is G — Go + GotGo- We assume that in tunneling 
the odd parity channel is fully transmitted (S 1 = —1) and the even parity one gives (see eq. (7)): 

t ~ VGd,dV, Gd,d = 



e - e Q + iT 



t° Iff-, §^=tan*° 

e - eo + iT ' Viet 



e- e 



. 2 tan 2 ^ T 2 

sin o = 5 — - = t ^t^ ?r- (17) 

l + tan 2 (5° (e-e ) 2 +r 2 V ; 

Here V is the tunneling matrix clement defined in eq. (2), eo is the rcnormalizcd quantity defined in cq. (11) and 
r = Tfi + Tl is the inverse lifetime of the resonance. This implies that eqs. (12,16) become: 

g = g u sm 5 =- , (18) 

(e- e ) 2 +T^ 

where g„ = (r 4r fp\ 2 . On the other hand, because StoG = — T/[(e — eo) 2 + T 2 ], eq. (4) becomes: 

where T = T^ri/r. At zero temperature both results of eqs. (18) and (19) coincide. 

At finite temperature, if the odd channel is fully transmitted cq.(12) can be generalized as 

g = g u J dcu ("^^) [-^(0)3m (t )] (20) 

where t is the t matrix above defined and related to the exact retarded Green function through the relation G = 
Go + GotGo- We have used the optical theorem: 

^ Tr {t^} = *le{f o t°} (21) 

which follows from the unitarity condition IS*! 2 = 1. 

III. ROLE OF THE ONSITE REPULSION U IN TUNNELING 

Up to now Ho referred just to a single impurity level e^. Indeed, charging energy U is the main feature of a QD 
and we have to deal with it. Something which is closer to a QD is an impurity level with onsite repulsion U. Inclusion 
of U in the Hamiltonian (1) leads to the single level Anderson model: 

Hand = -fficad + ed^^cr + t/^n CT n CT ' + y^ V a [c\ aa .da + d\c ka a]- (22) 

a- aa > kacr 

If Ed = —U/2 with respect to the chemical potential of the conduction electrons \x (taken as the zero of the single 
particle energies), the Anderson model which arises is symmetric. In fact, the energies of the empty impurity state, E, 
and of the doubly occupied impurity state, 2 E = 2ed + U, are both zero, while the singly occupied impurity level has 
energy 1 E = —U/2. 

In order to understand how the Coulomb repulsion on the dot site affects the Green's function of the dot we 
use a path integral formalism. We show that, in the limit in which the charge degree of freedom on the dot is 
frozen (U — > oo), the dot Green's function describes the dynamics of the only degree of freedom left, that is the dot 
magnetization < nj — 1/2 >=< 1/2 — nj > which is forced by a stochastic field X(t) in imaginary time, produced 
by the scattering of the lead electrons, with a gaussian probability distribution. The excitation energy associated to 
it is shifted from e^ to the Fermi level: this is the origin of the resonance at the Fermi level in the Kondo problem. 

As a first step, we will rephrase the previous result of an impurity with U = in this new approach. After linearizing 
the bands, the leads action can be written in terms of the field operators ip ua (where a = L, R and u includes all 
other quantum numbers) as: 



r0 r ^ 

S\ ea d = -iv F dr dx^2 X! ^ta( x i T )-i-' l Pua(x 1 T). (23) 

J J „, t n> 



" u a=L,R 



With respect to the tunneling action, if the L and R barriers are equal and the dot is zero dimensional the interaction 
term only includes the symmetric combinations 3>u(t) = —ts{4>ul{x = 0,t) + i/j u r(x = 0, r)) at the origin: 



r/3 
/o 
The total action is: 

/o 



S T = / dr>/2X)[V#t WduW + K*4(t)* u (t)]. (24) 

Jo 



^ = jf<ME<to 
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d u (r) } + Siead + S T . (25) 



One can first integrate out the degrees of freedom of the ip U a(r,x) fields for x ^= which are free-like, and next the 
field <&„ at x = which is the only one interacting with the dot variable d u (r). Because the fields in the leads are 
non interacting, the result of the gaussian integration yields an effective action for the dot: 

-S% s ex In J J J] l[v^ ua Z>VL e~ A | = -(3j2 d U^n) (iw„ - e d - E(iw„)) d„(iw n ) (26) 

I a=L,R u \ iuj n 

where the self-energy correction to the Green's function of the dot T,(iuj n ) = - l n ( ^"+g ) was obtained in eq. (10). 



A. Coulomb repulsion on the dot: freezing of the charge degree of freedom 

We now discuss the role of the onsite Coulomb interaction. In the large-t/ limit we have: 

f J r I i \ TT T -¥ f' 3 dT U n r+ n l) 2 + Tr i d(n r +n l )} 1L f d-r(n T -n, ) 2 

Jo 

= e^5(ni+n i +2e d /Uye%S« ^(™T-m) 2 } (27 ) 

where the delta function, which is defined by the last equality in the limit £/ — ► oo, implements the constraint of single 
site occupancy in the symmetric case, e<2 = —U/2. 

The quartic interaction is decoupled by means of a Hubbard-Stratonovitch boson field X(t), according to the 
identity: 

% J^drini-nif _ f T) X ~w /„" dr(X 2 (T)+2U(n- l -n l )X(r)) 



e ~Jo <"vn-»i; = / Me -«7 j 

Having introduced the field X(t), the partition function Z(fi) takes the form: 

Z( M )oc [ DXe~^ $o dTX2{T) 

x fm (DdiDdte-fo^'^G^-r')^ 

xS^+n.-l) 1 - J2 e { - 1)S ^ dT[n ^ ]X{T) . (28) 

j=T4 

Note that now the term ^dYlii n i was included in eq. (27), so in this case the dot Green's function is G7 Jioj n ) — 
iu> n — T,/Q\(iu> n ). At odds with the case U = here the resonance is at the Fermi level, in spite of the fact that 
the original localized level is at ed- The partition function in eq. (28) describes an effective spin-1/2 coupled to the 
fluctuating magnetic field X(t); its dynamics is constrained by the requirement that the impurity is singly occupied. 
It could be shown that the single site occupancy constraint is fulfilled in the average when the partition function of 
eq. (28) is used. 
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B. Quenching of the magnetic moment: singlet ground state 

The representation of the partition function given in eq. (28) allows us to recognize the doubly degenerate ground 
state (GS) of the impurity spin Sd = 1/2 with S^ = (— 1) CT (n CT — 1/2), a =f, J,, driven by the field X(t) and produced 
by virtual tunneling of electrons on and off the dot at energy /z. 

Anderson, Yuval and Hamann [21] calculated eq. (28) by integrating out the impurity (d,d) fields) and showing 
that the interaction with the delocalized electrons gives rise to a dynamics of the field X(t) which also interacts 
with itself at different times according to a logarithmic law. The problem was solved by mapping eq. (28) onto the 
partition function of a Coulomb gas (CG) [34] of flips in 1 — d, as we show in some detail in this subsection. 

Let us define the new field £(r) = X(t)/T, where T = 7w(0) \V\ is related to the tunneling strength. At low 
temperatures, saddle point solutions of the resulting single particle effective action are sequences of instantons 
£±( T >0 = ±CAftanh((r — t{)/tm) (where Ti,l — 1,2,... are the centers and tm ~ — is the width, that is some 
cut-off which regularizes the theory at short times) corresponding to jumps between the two minima of the effec- 
tive potential V e //[£] = (T 2 /ll) £ 2 - 2T/n [£tan _1 (£) - \ ln(l + £ 2 )] , located at £ M = ±U/2T, and interacting via a 
logarithmic potential a 2 In \(ri — tj)/tm\- 

For such a CG we can define the bare strength a 2 of the logarithmic interaction (the "charge") as a 2 = 2(1 — 8T/Uir) 
and the "fugacity" Y as Y = TMe~ Stot . 

Thus, the full partition function may be approximated with the sum over the trajectories given by hopping paths 
and will be written as: 

i-timf—- r" ^[^-.'-"""-""i^n (29) 

^ ( 2iV ) ! JO T M Jo T M 

that is the partition function of an effective one-dimensional gas of spin flips. The integral over the "centers of the 
instantons" has to be understood such that r, and Tj never become closer than tm- 

Now, we arc ready to perform the RG analysis to get the behavior of the model at large time scales (low temper- 
atures) [21]. The bare fugacity of the CG is Yj, = tm exp(— A) where A ~ tmU is the action of one single flip. The 
scaling of the fugacity and the renormalization of the coupling constant induced by processes of fusion of charges lead 
to the following Renormalization Group (RG) equations [34]: 

dY nt 2 drr 2 

= (l-^)Y;-^ = -2Y 2 a 2 (30) 
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(see Appendix A for a sketch on the derivation). 
We see clearly that the flow is towards Y — > oo and a 2 — ► and the scaling invariant energy (that is the Kondo 

temperature which we will introduce in the next section) is Tr- = t m e -1 ^ 1- ^"^ ~ (UT)^e^ 7TU ^ 8F \ 

Condensation of instantons in the doubly degenerate GS leads to the Kondo singlet < S z >= on the dot. Now 

we present an heuristic argument for such quenching of the total spin S on the dot; it runs as follows. 

The dynamics of the field £(r) with action A can be mimicked by a two- level system hamiltonian H21 with hopping 

energy A ~ \fn e ff [4) ~ \ 77- t m (^m/tm) 2 ~ U/2. The role of the interaction is to project out higher energy 

components from the dot states. Let us denote by |± > the two eigenstates of H21, i.e. the singlet and the triplet 
(with respect to the total spin of the dot plus conduction electrons, 5 = and S = 1) on the dot. Instantons, by 
flipping the impurity spin, produce a dynamics of the system "dot + conduction electrons " between these two states. 
Now, let us make an interesting analogy with thermodynamics. At temperature (3~ l the probabilities of having the 
system in one of the two states will be given by [35] : 

1 e- 2 ? x 

P+ = ^r : P- 



1 + e -2/3A ' " l + e -2/3A 

partly eliminating the component on the high energy eigenstatc | — >. Thus we are ready to give the connection with 
the CG picture described by the partition function of eq. (29). In our case the dynamics is given by the coherent 
zero point fluctuations of the system as a whole. From the statistical weight associated to a configuration with N 
instantons we can easily write down the formula: 

V°° 2ATK 2A 7f2An! d 

(N) = (N inst ) = ^"=0 OMI [ , = Y— ln(coshF) = ytanhF. (31) 

\ / \ mst/ ^ Y 2N (2N)\ dY y ' V ' 



The frequency A/2ir is obviously related to the average number of flips: f3\/2ir =< N >= TtanhY", as a direct 
implication of eq. (31). 

Hence, the probabilities of having the system in the states |± > within the ground state (T = 0) are: 

1 -i-rrYtanhY 

+ ~~ I _|_ e ~iTrY tanhY" ' = 1 _|_ g-47rytanhY ' ^ ' 

Because Y scales to infinity, the higher energy state completely decouples and the total spin on the dot is quenched: 
< S z >^ 0. 

IV. PERTURBATIVE ANALYSIS AT T » T K 

The Kondo Effect in metals containing magnetic impurities is responsible for the "anomalous" minimum in the 
resistivity p(T), as the temperature T ~ Tk' such a minimum in p(T) is due to scattering of conduction electrons off 
the localized magnetic impurities. Those contributions were first worked out by Kondo [36], who derived a correction 
Ap(T) ocln(^). 

As we pointed out in the introduction, a different realization of the Kondo effect may be achieved, in a controlled 
way, in a quantum dot in Coulomb Blockade (CB)-regime [3]. The Kondo effect in a dot is usually detected by 
measuring the linear conductance as a function of the gate voltage by connecting two electrodes to the dot. A dot at 
CB exhibits discrete energy levels. Changing the number of electrons is strongly prevented by electrostatic charging 
energy. Correspondingly, the total charge at the dot is quantized and the linear conductance is zero within large 
windows of variations of the gate voltage ("CB- valleys"). By analyzing the level structure of the dot it has been 
shown that, in some special cases, the dot at CB behaves as a magnetic moment antiferromagnetically interacting 
with the magnetic momenta of lead electrons. CB-valleys at different occupation number of the dot are bounded by 
resonant conduction peaks, as the chemical potential of the leads matches the one of the dot. Those peaks usually get 
sharper and better defined as T is lowered. However, as T <~ Tk, Kondo effect arises at the dot. Consequently, the 
CB-valley between two resonant conduction peaks "fills in" in a way that does not depend on the position of the Fermi 
level of the leads [3]. As T ~ Tk, the conductance g(T) exhibits the typical logarithmic raise [37]. The logarithmic 
raise, however, cannot hold all the way down to T = 0, as it must be limited by the unitarity bound. Therefore, a 
different approach is required in order to study Kondo effect in the T = 0-unitarity limit, as we will discuss in the 
next section. 

In this Section we derive the Kondo Hamiltonian from the single impurity Anderson model (eq. 22) focusing for 
simplicity on the isotropic case Jj_ — J z — J, then we sketch the pcrturbative Renormalization Group (RG) flow for 
the coupling strength of such a model. In general, the lower is T, the more spin-flip processes become likely, which 
make the running coupling constant J grow. At the Kondo temperature Tk the perturbative analysis breaks down, 
that is v (0) J{Tk) ~ 1, so that Tk is the characteristic scale that divides the regions of weak and strong coupling. 

A. Derivation of the spin dynamics hamiltonian 

In the following we will restrict our analysis to a two-fold degenerate dot level: in such a case the QD can be 
modelizcd as a spin-1/2 magnetic impurity. Indeed, because the charge degree of freedom is frozen at Coulomb 
Blockade we mimick the dot with its total spin Sd and describe its interaction with the delocalized electrons by means 
of the Kondo Hamiltonian: 

Htff = #lcad + Hk = -fflcad — JSdS(Q), (33) 

where <?(0) = X^fcfc' S ct <t' ( c I<t^ctct' c fc' V ) pl & y s the role of a spin density of the itinerant electrons at the impurity site 
(/? ' a<J i are the Pauli spin— | matrices). Its components behave as a spin— i provided single occupancy of the site 
x = is guaranteed. 

In the anisotropic Kondo model the couplings of the x, y components J±_ are different from the one of the z 
component J z . This effective Hamiltonian is more suitable to describe the low-T physics of the system because 
its dynamics shows how the system flows towards the nonperturbative regime. In the following sections by "low 
temperature" we'll mean T ~ Tk- 



As a first step we show that it is possible to derive the effective Hamiltonian (33) from the one for the single 
impurity Anderson model (eq. 22), where H& — £d^ a n a + UJ2 aa ' n a n a ' > by means of the Schrieffer- Wolff (SW) 
transformation [38] . More details can be found on the book by Hewson [5] . 

Let S be the space spanned by the (almost) degenerate dot's states. As the number of electrons at the QD is fixed, 
the relevant degrees of freedom can be described by an effective Hamiltonian H e g acting onto S only. In order to 
construct H c g, we apply the SW transformation to the Hamiltonian in eq. (22). We denote by P the projector onto 
S and by ed the energy of the states within 5, so that the effective Kondo Hamiltonian is given by: 



SH cS » P V(l - P) 



1 



£d — Hy> — H] 



lead 



-(1-P) V P. 



(34) 



Indeed, by inserting eq. (22) into eq. (34), we get the result: 

#cff = #lead - V (0) J2 — rfUr + E V » 

aa aa:k,k 



E v a v a , 

ot.a. \<J.a :k,k 



l 



i 



£d + U e d 



£d + U e d 



acr^ d ' T aa l C k' 



kao k aa 



(35) 



are the impurity spin components. Spin commutation relations 



where S\ = ^ CT aS a d a , S'j" = dldi and S^ = d]c 
are obtained if no double occupancy is admitted. 

Besides i^icad, the first and the second term at the r.h.s. of eq. (35) represent, respectively, a renormalization 
of the dot's energies and a potential scattering term; the third term is the one which induces spin flips. The two 
contributions appearing in the potential scattering term and in the spin-flip term refer to two virtual particle and hole 
processes, respectively. In the first case a particle is added to the dot so that the energy ed + U is involved, while in 
the second case a particle is subtracted from the dot level thus paying the energy |ed|. The potential scattering term 
vanishes if the Anderson model is symmetric, as we have assumed up to now [td = —17/2). 



B. Perturbative Renormalization Group Approach 



In this Subsection the starting point of our reasoning is, for simplicity, the isotropic limit of Kondo hamiltonian in 
eq. (33). Scattering of conduction electrons by the impurity produces a self energy correction, as well as a correction 
to the interaction vertex. The transitions between the states close to the Fermi level tF and the states within a 
narrow strip of energies of width SD near the edges of the band of width 2D are associated with an high energy 
deficit. Such transitions are virtual and their influence on the states near ep can be taken into account perturbatively 
in the second order. In the poor man's scaling approach one includes second order corrections arising from processes 
in which the electrons k! are scattered to an intermediate state at energy D > e&» > D — SD or — D < tk" < —D + SD 
where D is some ultraviolet cut-off. Because they involve high intermediate energies, one can think of including them 
perturbatively, so for an electron k' scattered into the final state k we have: 



E <k\H K \k" 

fe"G7 



>< k" 



1 



E — i^lcad 



w 



>< k"\H K \k' >~ -jvo)«^&-EE c L? 

kk era 



'C k ' 



(36) 



where 7 is the k" domain mentioned above. To justify the second step one observes that, in the shell 7, the operator 
Hicad can be replaced by an energy D and, in comparison to it, the eigenvalue E can be put at the Fermi level E — 0. 
The result is an effective Hamiltonian acting within the band of a reduced width D — SD, of the same form as Hk in 
eq. (33) but with a modified value of J . This gives the following correction to the previous coupling constant in Hk'- 
SJ ~ —J 2 v(Q)8D/D [39]. Successive reductions of the bandwidth by small steps SD can be viewed as a continuous 
process during which the initial Hamiltonian is trasformed to an effective low-energy Hamiltonian acting within the 
band of reduced width D — SD. The evolution of the exchange amplitude J during such a poor man's scaling can be 
cast into the form of a flow differential equation [39] : 



dJ 



-v(0)J 2 



d\nD 
Integration of eq. (37) gives rise to the renormalization of J: 



(37) 
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m-jh= ln w (38) 

where j = Jv(0), that is 

Scaling can be performed down to D ~ fc^T, also we see that, if j'q = i(-Do) > (antifcrromagnetic coupling), the 
running coupling constant j increases. Eq. (38) shows that De^> is a constant of this flow, which defines an energy 
scale as Tk'- 

k B T K = D Q e-^. (40) 

In general, at the Kondo temperature Tk the system has approached the scale at which the perturbative analysis 
breaks down, so that for T < Tk j starts to diverge in the flow. 

Now we discuss the conductance in such a perturbative limit. Let us take a S(x) -like dot with an even barrier 
potential (Tr = Tj,). To lowest perturbative order the leading self-energy correction to t°(u> + i0 + ) (see eqs. (16), 
(17)) (lu = 0) yields: 



!<() 



Vte{^t°\~ir 2 (v(0)J) 2 , (41) 

Kq 



where -77777 = In ■£- is the invariant charge defined through eq. (40), so that the conductance can be written as: 

g « 7T 2 In" 2 £-, T K <^T<^ D. (42) 

The tail of the conductance in the perturbative limit T 3> Tk is logarithmic. 

C. One channel anisotropic Kondo model and the Toulouse limit 

In the following we rewrite the Coulomb gas approach introduced in subsection V.B focusing on the single channel 
anisotropic Kondo system. Let us start from the general effective Hamiltonian: 

F e A ff = H lcad + H K = tficad + JzSfaiO) + ^ (S+cj- (0) + Sja+(0)) . (43) 

As stated in ref. [21], the perturbation term (the one containing J±) has the effect of flipping the local spin at 
each application; hence the problem of calculating the partition function of such a system reduces to the evaluation 
of the amplitude for a succession of spin flips at times ti,T2, ... and the Feynman sum over histories is the sum over 
all possible numbers and positions of flips. So, we get the expression (29) which can be rewritten as: 



^,r ,2N f dT2N f T2 -™ dn +£ 2 " (_l)«+i(2- e) l„|^i| .... 

jcg = > (J±T M ) / •■• / [e ^»^=i I tm I]. (44) 

JO T M Jo T M 



N=0 



Here e is the quantity: 



8S AF 8S\ F ^ 

5— — UzTM (.45) 



where Saf is the scattering phase shift of antiferromagnetic sign introduced by the J z 5Jct z (0) term. It can be seen 
clearly that the sign of J z , or e, determines whether the coupling is ferromagnetic or antiferromagnetic, so the condition 
e > gives the antiferromagnetic coupling. 

Let us now notice that eq. (44) is a function of three parameters only: -^— , Jj_tm,£', we are interested to the case 
t — > 00 (low temperature). In such a case it is well known that the ferromagnetic Kondo system has a mean spin 
moment which corresponds to a long range order of the classical system: the positive and negative charges are all 
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bound in pairs forming dipoles all pointing in the same direction. The phase transition occurs, at least for small J±, 
at the ferromagnetic-antifcrromagnctic boundary point e = 0. 

For such a system it is possible to derive a set of scaling equations by renormalization of the cut-off tm] these laws 
are exact for small Jj_tm and e. The physical picture is that of many close pairs of flips which change slightly the mean 
magnetization, located between pairs of isolated flips which are reversals of the magnetization over a larger timescale; 
so, the isolated flips can be considered as acting in a medium where the close pairs modify the mean magnetization. 
This line of reasoning leads to the following scaling laws for the "fugacity" J±tm and the "charge" e respectively: 

d -^ 1 = \ GW) ; -£- = (2 - e) [J,r M f . (46) 

a In tm 2 dm tm 

Let us make a few comments on such equations. First, these are also exact for finite e, because only Jj_tm need 
to be small. Second, they are compatible in the isotropic case J z = Jj_, s ~ 2J z tm — 2J±tm where J±tm and e are 
small; because J \Tm and e scale together, the isotropic Kondo system remains isotropic at every time scale. In the 
anisotropic model, eqs. (46) become: 

</:- . ^J^TM . £ 2 _ 4J 2 t 2 m = mnst ^ (4?) 



d(J±T M ) 



thus the scaling lines are a set of hyperbolas with asymptotes corresponding to the isotropic case. 

We see clearly from Fig. 3 that all ferromagnetic cases below the isotropic one scale onto the case J\_tm — which 
is an ordered one, so the transition line coincide with the ferromagnetic case. Conversely, in the antiferromagnetic 
(AF) case J±tm and e increase starting from some values, small at will, and we always can find a timescale for 
which £ ~ 1. Such a timescale is a crucial one for the Kondo phenomenon, it is the Kondo temperature already 
introduced in the previous sections and sets the scale at which the system behaves as it were strongly coupled. Thus, 
the renormalization procedure is valid up to e ~ 1. 

Now, let us briefly focus on the so-called Toulouse limit J z = J±, s ~ 1; in such a case the system is equivalent to 
a simple quadratic model with Hamiltonian: 

H T =J2^c{c k + VY^[^c k + cld}. (48) 

k k 

Here Ck,c k are Fermi operators for free spinless electrons and d,dft are Fermi operators for a local resonant state. 
The partition function corresponding to the Hamiltonian (48) is exactly equal to the one in eq. (44) for the case of 
classical charges ±1. Such a theory corresponds to free particles and there is no renormalization. 

V. BREAKDOWN OF PERTURB ATIVE APPROACH ATT~T K : STRONG-COUPLING FIXED POINT 

In the previous Section we have derived the perturbative flow equation for the Kondo coupling constant J . Integra- 
tion of the renormalization group (RG) equation shows that J grows as T is lowered, until the perturbative analysis 
does not make any sense anymore. The question whether the RG flow stops at some finite-T scale fixed coupling 
or goes all the way down to T = has been widely debated in the literature (see [5] for a review on the subject). 
However, from numerical RG analysis and from the exact Bethe-ansatz solution of the model [7] it is now clear that 
the system takes no fixed points at finite T, but the RG flow goes all the way down to T = 0. So, the aim of this 
Section is to discuss the physics of the Kondo system in the T = 0-unitarity limit. 

In the case of a localized spin- 1/2 impurity antiferromagnetically interacting with the spin of one type of itinerant 
electrons only (spin-1/2, one channel Kondo effect), the theory of the unitary limit was first developed by Nozieres 
[40]. As T is lowered all the way down to 0, the flow of the coupling strength runs all the way toward an oo-coupling 
fixed point. At the fixed point, the impurity spin is fully screened and the localized magnetic moment is effectively 
substituted by a spinless potential scattering center with infinite strength (which forbids double occupancy at the 
impurity site). The T — > "unitarity limit" is well-described by a Fermi liquid theory, which also allows for calculation 
of finitc-T corrections to fixed-point values of the physical quantities. Both elastic and inelastic scattering processes 
provide finite-T contributions to the conductance. As we will show in the next section, the two kinds of processes 

contribute at the same order, thus giving raise to corrections to the unitary limit proportional to I -i^— ) [40]. 

A generalization of Kondo's original idea was done by Nozieres and Blandin [41]. They showed that interaction of 
itinerant electrons with magnetic impurities in metals may involve electrons with different quantum numbers besides 
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the spin: for instance, electrons with different orbital angular momentum. Different orbital quantum numbers define 
different "channels" of interaction. Therefore, "many channel'' Kondo effect may arise. Although in the perturbative 
region there is no qualitative difference between one-channel and many-channel effect, deep differences may arise in 
the unitarity limit, depending on the number of channels K and on the total spin of the impurity, 5. 

In the case 25 = K, at T = the magnetic moment at the impurity is fully screened by electrons from the leads. 
The impurity forms a singlet state with 25 conduction electrons and no other electrons can access the impurity 
site. The system behaves exactly as if there were no impurity, besides a boundary condition on the wavefunction 
of conduction electrons which takes into account that the impurity site is "forbidden" to them [40]. Such a Fermi 
liquid state is stable against leading finite-T corrections, as we shall see below, and corresponds to an infinite effective 
coupling v (0) J. A special case of this is the one channel Kondo just discussed, with 5 = \ and K = 1. 

In the case 25 > K, in the strong coupling regime a residual magnetic moment is still present at the impurity, 
since there are no more conduction electrons able to screen the localized spin. The corresponding fixed point is again 
a Fermi liquid, but with a localized partially screened magnetic moment at the impurity that is ferromagnetically 
coupled to the itinerant electrons: it again provides the stability of the local Fermi liquid. 

A very special case is the 25 < K one. In such a case conduction electrons attempt to "ovcrscrccn" the impurity 
in the strong coupling limit, that is the resulting magnetic moment is opposite to the original one. The coupling 
among the localized residual magnetic moment and the itinerant electrons is now antiferromagnetic. It drives the 
system out of the strongly coupled regime towards a finite-coupling fixed point J* [41], as we clearly see in Fig. 4. At 
J = J* the leading finite-T interaction is given by a marginal three-particle operator which breaks down the Fermi 
liquid state and generates a non- Fermi liquid behavior in the physical quantities. We will discuss such an issue in the 
Section VII of this review. 

VI. THE FULLY SCREENED SINGLE CHANNEL KONDO CASE 
A. Finite temperature corrections to the one-channel Kondo conductance 

At T = the linear response of the Kondo system to an applied voltage bias reaches the so called unitarity limit. 

The response function is the resistivity in a bulk system (magnetic impurities in diluted alloys), while it is the 
conductance in a quasi one-dimensional system as the system of our interest: a dot with applied contacts. 

The striking feature of the Kondo effect in diluted alloys is the minimum in the resistivity at low temperature which 
violates the expected property: dp/ dT > 0. Indeed, well below Tk the resistivity increases again up to a maximum 
value proportional to the number of impurities Ni per unit volume ("unitarity limit"). 

On the contrary a maximum of the conductance is expected at T = for Kondo conductance across a quantum 

Op 2 

dot, where the unitarity limit that is reached is —g u - 

Being the s— wave scattering effectively one-dimensional like, such a 'reversed behavior' looks paradoxical. However, 
it is just a consequence of the difference between the 3-d and 1-d impurity scattering, as we explain here. In order to 
illustrate the difference, it is enough to note that both facts stem from the main feature of Kondo impurity scattering: 
the formation of a resonance at the Fermi level due to many-body effects, which implies that the phase shift reaches 
the value S = tt/2 and that scattering is resonant at the impurity. This produces different results in the two cases: 

3-d: a spherical s— wave is diffracted from the impurity with maximum amplitude at the resonance, what increases 
the flux propagating backward and enhances resistivity. 

1-d : resonant scattering coincides with resonant trasmission in this case, what implies the vanishing of backward 
reflection and enhances the conductance. 

Temperature corrections are twofold. One is due to the energy dependence of the phase shift close to the resonance 
S°(e) = tt/2 — ae 2 and to the fact that energies close to the Fermi surface are sampled, because the Fermi functions 
are not step-like at finite T. The second one is due to inelastic processes which produce transitions from the singlet 
ground state to the excited states. The latter can be accounted for with an expansion in inverse powers of the singlet 
binding energy [40] . We include a simplified approach to the problem which rests on the Fermi liquid nature of the 
excitation spectrum in Appendix B [40] [8] [10]. It is found that corrections to the conductance are quadratic in 
temperature, as it is usual in the Fermi liquid theory: 



9 = 9u 






T<T K . (49) 
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The weak-coupling (T >• Tk) and the strong-coupling (T <C Tk) asymptotes of the conductance, eqs. (42) and (16), 
have a very different structure but, since the Kondo effect is a crossover phenomenon rather than a phase transition 
[6] [7], the dependence g (T) is a smooth function [10] throughout the crossover region T ~ Tk'- 

9 = 9uf (J^ ■ (50) 

The universal function f(x), as found by resorting to numerical renormalization group (NRG) in refs. [10] [11], is 
plotted in Fig. 5. It interpolates between / [x ^> 1) = ^g- (In a;) and / (0) = 1. 

B. The Kondo resonance in a magnetic field 

A small magnetic field B lifts the degeneracy of the spin states at the impurity. This produces a splitting of the 
resonance and a change of its shape which has been numerically studied mostly via NRG [11]. The splitting of the 
peak increases linearly with B and is A = 2ugB, i.e. twice the Zeeman spin splitting. 

This can be understood easily if one considers the particle and hole virtual occupations which mediate the Kondo 
interaction. 

Let us consider the symmetric case ed = —U/2. The dot states with N electrons and one unpaired spin are N | 
and N J. corresponding to e^ + B/2 (u B = 1) and e^ — B/2, respectively. In presence of spin splitting two particle 
processes are allowed: p a {p-a) in which a spin — a is added to the state Na and subsequently a spin — a is removed 
with the energy balance: 

Pi -» 5E(N t, + 1) + 8E(N + i,-l) = U + (-U-B) = -B 

pj -» SE(N |, + T) + SE{N + 1, - |) = (U + B) + (-U) = B. 

Similarly the hole processes have an energy balance: 

fr T - SE(N t, - T) + SE(N - 1, + 1) = {-e d - B/2) + (e d - B/2) = -B 

h i -» 5E(N |,-|) + SE(N - 1, + t) = (-e d + B/2) + (e d + B/2) = B. 

This implies that there are two peaks in the spectral density at w = ±B corresponding to the pi, hi and p^, hi spin 
flip processes, respectively. 

The value at the Fermi energy of the spectral function is related to the B-depcndcnt phase shift in the Fermi liquid 
picture by the Friedel sum rule: 

-txv (0) 3?m {t a (u) = 0, T = 0, B)) = sin 2 5 a (B). (51) 

The Bethe Ansatz solution of the problem relates the phase shift at the Fermi level to the magnetization of the 
impurity: S a (B) = f [1 — 2Md(B)] [43]. Hence a reduction of the peak height with increasing of the magnetic field 
follows. 

It has also been argued that the renormalizability of the Kondo problem could break down when B exceeds some 
critical value related to the Kondo temperature, because of the back action of the induced conduction electron 
polarization cloud on the impurity leading to a broken symmetry state with < S z >^ [44]. 

C. Crossing of the dot levels in a magnetic field and enhancement of the Kondo conductance 

Conventional Kondo resonant transmission requires a magnetic moment to be present on the dot. Usually dots 
with an even number of electrons are in a singlet state, while dots with an odd number of electrons have an unpaired 
spin and have a doublet GS. Hence, there is a parity effect: CB conduction valleys with N = even do not display 
the Kondo conductance while those with N = odd do. An exception to this rule occurs at zero magnetic field when 
Hund's rule applies. This was found experimentally in vertical QD e.g. at TV = 6. The GS of the isolated dot has 
S = 1 (triplet) [45]. An underscreened Kondo effect is expected at this point and the GS of the interacting system 
becomes a doublet. By applying a weak magnetic field orthogonal to the dot (B = Q.22Tesla), a transition of the GS 
from triplet to singlet (T- S) has been found. The single particle energy levels become angular momentum dependent 
and Hund's rule no longer applies. Close to the crossing a remarkable enhancement of Kondo coupling with increase 
of the Kondo temperature was experimentally found. Indeed, scaling shows a non universal critical temperature when 
the interplay of the fourfold degenerate states (\SS Z > with S = 0, 1) is included [46]. 

An extended and unified approach to the problem can be found in [47]. A minor difference is the fact that they 
consider an in-plane magnetic field as the source of the crossing which could produce Zeeman spin splitting (ZSS) of 
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the single particle states. The four dot states are mapped onto a two impurities Kondo model (2IKM) [48] with spin 
Si,D2- However they are coupled by a potential term Vnp nn (0) £> r S 2 and an exchange term ilnsnn{0)-o iX~S 2- 
Here p, a are the charge and spin density of the conduction electrons with n and n 7^ n labeling two different orbital 
states. Because these terms violate the invariance under particle-hole transformation, the 2IKM cannot flow by scaling 
to the non Fermi-liquid fixed point, which is known to be a remarkable feature of the model. In the case of large ZSS 
A the RG flow terminates at D w A. Two of the four states are ruled out in the flow and conduction electrons couple 
to one single effective spin 1/2 with one extra unusual term in the effective Hamiltonian which is a Zeeman term for 
the conduction electrons. 

Kondo conductance can take place in a dot with an even number of electrons also in a strong vertical magnetic 
field [35]. Orbital effects induced by B can produce the reversed transition from the singlet to the triplet state (S-T). 
Indeed, a vertical magnetic field on an isolated dot favors transitions to higher spin states [49]. In this case the ZSS 
is anyhow sizeable and the crossing involves the singlet state and the component of the triplet state lowest lying in 
energy. In a vertical geometry with cylindrical symmetry orbital angular momentum to and z— spin component a are 
good quantum numbers. In particular because the singlet state has total angular momentum M = and the triplet 
state involved in the crossing has M = 1, only a (m = J.) electron can enter the dot when it is in the triplet state. 
On the contrary, only an (to = 1 f ) electron can enter the dot, if it is in the singlet state. This implies that there is 
one single channel of conduction electrons involved in cotunncling processes, with orbital and spin degrees of freedom 
locked together. Again, a residual effective spin 1/2 survives at the dot, while N is even [23]. This is another striking 
manifestation of the spin-charge separation that occurs at the Kondo fixed point. Usual Kondo coupling leads to 
N = odd together with 5 = 0. In the complementary situation here described, it is N = even and S = 1/2. 

VII. THE OVERSCREENED TWO CHANNEL KONDO CASE 

A two channel Kondo behavior has been invoked in an experiment by Ralph and Buhrman [50] on clean Cu point 
contacts, defects in the metal that can be described by two level systems (TLS). The TLS could tunnel between 
the even and odd state with the assistance of the conduction electrons. Their physical spin is not involved in the 
scattering, so that two channels arc available [51]. It is still unclear whether the Kondo temperature can be large 
enough so that any effect of the Kondo correlation can be measured [52]. However, these experiments have triggered 
renewed interest in two channel Kondo conductance. The temperature and voltage dependence of the conductivity 
have been numerically calculated within the "Non Crossing Approximation" (NCA) [53] and found to be consistent 
with a scaling Ansatz motivated by the Conformal Field Theory (CFT) solution of the problem [8]. Since then, no 
other experimental proof of two channel Kondo effect in impurities has been produced. The result for the imaginary 
part of the transmission t is: 

9m*(w) - Sjf. (52) 



From this result, we get a \JtJTk temperature dependence of the conductance as T —> 0, which is a clear signature 
of the breakdown of the Fermi liquid [8] (see Fig. 6). 

In order to emphasize the deep difference between single-channel and many-channel Kondo effect in the T = 
limit we just mention here that the imaginary part of the proper self-energy, close to the Fermi surface, behaves 
as SmE(fc,o;) ex C^lo 2 (where the chemical potential a = is taken as reference energy) in the Fermi liquid case 
and in the single channel Kondo problem. On the contrary, it behaves as $JmS(fc, u>) oc C" k u>z in the two channel 
" overscreened" Kondo problem. In the following we refer, for simplicity, to the two channel " overscreened" case. 
Through the Kramcrs-Kronig relation 



d 1 r 00 
— KeE(fc,u;) = P / duj 

duj TT J_ QO 



,d u ?smT,(k,u)') 



we see that ^j5fteE(fc, w)| w ^o is finite in the first case, at the Fermi surface, on the contrary it has a power law 
divergence in the second case. It follows that the quasiparticle pole residue Zk = [1 — 9wE(fc,w)] | w _>o vanishes as a 
power law in the second case at the Fermi surface. Luttinger [54] showed that, to all orders of perturbation theory 
in the interaction, the imaginary part of the proper self-energy behaves as $JmE(fc, uj) oc C^uj 2 {Cf. > 0) close to the 
Fermi surface, which implies an infinite lifetime for the quasiparticles at the Fermi surface. The Fermi surface is sharp 
and well defined. These are the foundation stones of the normal Fermi liquid Theory and they are invalidated in the 
spin-i two channel Kondo case. 
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In this Section we focus, in particular, on two channel spin-1/2 Kondo effect in both the perturbative region and the 
unitarity limit. We describe the scaling perturbative approach and the bosonization (T ~ 0) technique, respectively. 
An attempt to extend the Anderson- Yuval approach of Section IV. C to the two channel Kondo case can be found in 

[42]. 

A. Perturbative analysis at T >> Tk 



The starting point of our reasoning is the effective Hamiltonian given in cq. (35), which we will take in the isotropic 
limit. In the following we will restrict our analysis to a two-fold degenerate dot level. In this case, as we pointed out, 
the QD can be modelized as a spin-1/2 magnetic impurity, whose spin is given by: 

S a d = ld!,T a ,d> 

« 2 7 77 7 

(that Sd is a spin-1/2 comes out from the identity S% = 3/4, valid in the case of single-occupancy for the dot's level). 
In this case and for a generic number of channels for the itinerant electrons, eq. (35) takes the form: 

^ = J E (e«1^V)E E (cL^ava)) ( 54 ) 

a \77 kk era ;a / 

where a G (1, ..,K) is the channel index and the constant J is taken as a perturbative parameter (> 0). Infrared 
divergent diagrams provide a flow of J as a function of T. We are now going to derive the perturbative /3-function at 
third order in J . At finite T the Green function in Fourier space will depend on the momentum of the particles and 
on the Matsubara frequencies u> m — M-{m + i) (for fermions). In our case, the Green function for the lead electrons 
is given by: 

CO. „„' (*w mi k) = FT \(f[c aa (r, k )cj, ,(0, *)]>} = AA- (55) 

I a a ) iu) m — vpk 



aa ;aa 



where FT stands for 'Fourier transform' and T is the time ordering operator, while the Green function for the 
<i-fcrmion is: 

G W (iu> m ) = FT{(f[d 7 (r)4(0)})\ = 5 -^L. (56) 

The interaction vertex determined by H e e is: 

V ao i' .{{ilo^}) = 6<-T a ,r a ,tf( w W +w (2) _ w (3) _ W W). (57) 

(T(T '77 ^ ^ *^ -* ' QOL A 7^y (7(7 V Tfl Tfl VTl ITL / \ I 

The one-loop structure of the theory provides a renormalization to the interaction vertex, that is, to the coupling 
constant, by means of the two diagrams drawn in Fig. 7. The sums over Matsubara frequencies can be calculated by 
using the standard techniques described, for example, in [55]. In the low-energy limit for the electrons from the leads 
(that is, if only excitations about the Fermi level are taken into account), the sum of the two diagrams is given by: 

©i +V2 * -i^{[8 aa ,5 a0i> +T" aa ,T a aa ,} + [-<U«U< +r CT Vr Q V]}2,(0)ln(^) (58) 

(D is an ultraviolet cutoff, identified with the width of the conduction band). The corresponding renormalization to 
the coupling constant J is easily worked out and is given by: 

A( 2 )j(T,^) = 2^0)J 2 ln(-^-). (59) 

k B T 

A careful analysis of the vertex renormalization at third order in J reveals that several third-order diagrams are 
already taken into account by the scaling equation generated by the second-order vertex correction, as discussed in 
[56]. The only "new" contribution comes from the "non-parquet" diagram shown in Fig. 8. Because of the loop over 
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the fermions from the contacts, such a contribution carries an overall factor of K, that is, the diagram is proportional 
to the number of channels and the correction to the coupling constant up to third-order in J is: 

A< 3 > J = (21/ (0) J 2 2K (u (0)) 2 J 3 ) ln(-r^=). (60) 



Integration of the renormalization group equation for J provides: 

2. 7 2j 2 n [j \l~Kj 



= x - x (61) 



where j = v (0) J and x = ln(^^). 

It is seen from eq. (61) that the fixed point is reached at an intermediate coupling j* . However, the perturbative RG 
analysis is not conclusive here, because even in the ordinary, single channel, Kondo model an artificial intermediate 
coupling fixed point is produced when the perturbative RG equations are expanded to 3rd order [56]. 

It is straightforward to infer the Kondo temperature Tk from cq. (61). Usually Tk is defined as the temperature 
scale at which j becomes 0(1). Following such a criterion, at T = Tk we can neglect 1/j compared to 1/jo and 
obtain the following approximate formula for Tk'- 

k B T K ^D(j )^e-^. (62) 

Eq. (62) is quite general, in that it provides the value for the crossover temperature for any number of channels 
K. This proves that the way the system approaches the scale at which the perturbative analysis breaks down does 
not depend on the number of channels, except for a redefinition of the Kondo temperature. Hence, logarithmic 
divergencies are expected when approaching to the Tk, no matter on what the number of channels is. In the next 
subsection we will analyze the T = behavior of such a system, and will see that it is dramatically different from 
the one channel case for what concerns the fixed-point properties. 

B. Analysis at T ~ 

Several techniques have been applied in order to get informations about physical quantities around the fixed point 
in the Kondo regime in such a limit. Finite- T corrections have been derived by means of bosonization techniques 
[57] [58], of Bethe-ansatz like exact solutions [7] and of Conformal Field Theory (CFT) techniques [8]. Now, CFT 
approach is extremely effective in calculating finitc-T corrections, Wilson ratios and other exact results concerning 
Green functions [8], but its starting point that charge, spin and "flavor" quantum numbers are enough to identify a 
primary fcrmionic field in the theory leads to an inconsistency: the corresponding on-shell S matrix comes out to be 
0. The solution to such an "unitarity paradox" has been suggested by Ludwig and Maldacena [59], who introduced 
a fourth " spin- flavor" quantum number. So, while at the unitarity limit charge, spin and flavor do not change upon 
scattering off the impurity, the spin-flavor changes, giving rise to one more scattering channel. Such an extra quantum 
number allows for an off-diagonal on-shell S matrix, therefore the unitarity paradox simply means that the diagonal 
elements of the S matrix with respect to the spin-flavor number are 0. However, a theory of the unitarity limit, needed 
in order to compute, for instance, within an unified framework transport properties with finite-T corrections, is not 
yet fully developed. In this subsection we briefly sketch the first steps in order to derive an appropriate scattering 
potential a' la Nozieres, in the case of two channel spin-1/2 overscreened Kondo effect. We follow an approach 
equivalent to the one by Tsvelick [60], that is the introduction of a rcgularization procedure able to move the fixed 
point toward infinite-coupling. Then, we go through a bosonization - refcrmionization in order to account for the 
scattering processes with changing of the spin-flavor. The bosonization procedure allows us to split the degrees of 
freedom involved in our problem: in this way we will show that the Kondo interaction involves only the spin and 
spin-flavor degrees of freedom while the charge and flavor ones are fully decoupled. This remark is a crucial one 
because it makes possible to derive a S matrix in the unitarity limit which comes out to be diagonal in any quantum 
number, but the spin-flavor; such a S matrix describes the whole system of dot and contacts. 

In the following we discuss the case K = 2, S = 1/2, so we have two flavors of conduction electrons from an 
effectively one-dimensional conductor which interact with a localized spin— 1/2 magnetic impurity. 

Let c aa (x) be the lead electron operators (o —] , { is the spin index, a = 1, 2 is the flavor index). A lattice- model 
complete Hamiltonian is written as: 
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H 2Ch = -t^y^i&Ccvix + a) + c) aa {x + a)cw(x)] 

X 

-^J2 c ^( x )c a Ax) + J^ d -^i(0) + ^m (63) 

X 

where a a (x) = \c i a ak x )^aa' c aa' ( x ) an< i a i s the lattice spacing. 

Now, we can linearize the dispersion relation around the Fermi surface and introduce two chiral fields c±_ aa {x). 
Even and odd parities can be introduced to obtain fields with the same chirality, (\f aa and 4>° aa , so odd parity fully 
decouples from the interaction Hamiltonian. In order to properly deal with the interacting fields, we will bosonize <\> e aa \ 
in particular, we define four bosonic fields f M in terms of which it is possible to express densities of the remarkable 
physical quantities. To this end we are led to construct the four bosonic fields ^>x (X =ch,sp,fl,sf) for the charge, 
spin, flavor and spin-flavor degrees of freedom [58] [59] [61] as linear combinations of the previous ones. 

In this way it is possible to realize two "inequivalent" representations of the fields <f> e aa in terms of the fields ^ x 
(X =ch,sp,fi,sf), (j) e aa .j and C CT;/7 , given by: 

<t>la;l( X ) = V™ ■ e-5[*ch+^ sp +«* fl +a CT * af ](*) . 

where r\ and £ are suitable Klein factors. Notice that the two fields differ only in the spin- flavor quantum number but 
such a difference is a crucial one. 

Now we are ready to rewrite the two channel Kondo interaction Hamiltonian in bosonic coordinates ^x as: 



H 2 K ch = jU+: e -^p(°) :: cos(tf sf (0)) : +SJ : e 4 *- (0) :: cos(* sf (0)) : +5|- 



L dtf sp (0) 



'2-k dx 
J^ d -[S g (0) + S„(0)] (65) 



where the spin densities IU,b(x) are given by: 



Ait dx ' 



E±(x) = : e^^-p+^Kx) . Sff (x) = __[* Bp + * Bt ](x), 



£±(x) = : e ±*[*.p-*-](-) : £*(*) = ~[* sp - *■£](*)■ (66) 

Both £ 9 (x) and S u (x) are SU(2) spin current operators, so the vector space is made of the bosonic vacuum \bvac) 
and the bosonic spin— i spinors at a point x: 

\<j g ) = : e^^p+^K*) : \bvac) 

\a u ) = : e ^[* a p-* S f](*) . \ bvac y ( 67 ) 

No triplet adding up g and u spin species together can occur because, given S ff)U = J_ L / 2 dx : £ ff)U (a;), we have: 

E s K)=0 S„|c7 g )=0, (68) 

that is, if at a point x the spin density associated to S s is ^ 0, then, at the same point, the spin density associated 
to S M is = 0, and vice versa; this does not allow for overscrccning. Such a statement is a crucial one, indeed it 
corresponds to a particular rcgularization scheme [60] able to move the finite-coupling fixed point corresponding to 
the unitarity limit down to J = +oo. At this infinitcly-strongly coupled fixed point the impurity spin will be fully 
screened in a localized spin singlet. In principle, the system might lay within any linear combination mixing the two 
representations <?, u, of the form: 



|GS)„| x = =|ft>®^( !,/-/<! I„) 



-|*>®£(IT fl > + /'l D 

x=0 Z 



(69) 

x=0 



where \ it), | -IJ-) are the two impurity states with opposite spin polarizations. Then we search for the two independent 
linear combinations that do not change upon scattering of lead electrons; it can be shown that such combinations 
correspond to the values jj, = ±i. 



Now, in the fixed point limit, we " refermionize" . Physical states mix both representations, I and II. Scattering 
by the impurity states should conserve all physical quantum numbers. It can be shown that the elastic scattering 
in the unitarity limit swaps the two inequivalent representations I — II of the lead electrons. Correspondingly, the 
impurity absorbs/emits one spin-flavor quantum. The even parity S matrix has the following representation in the 
I — II space: 

sv=o)=(° -;).( e ? f *-;;). <™> 

According to eq. (70) the phase shifts induced by the scattering are S 1 II {uo = 0) = — f , 8 n ! (oj = 0) = f, while in 
the case of one channel spin— 1/2 Kondo effect the phase shift was 6° = -|. Finally, the conductance can be easily 
obtained by the Landauer formula using eq. (70) and S 1 = — 1: 



g = Tr {T} = Tr 



^ 




(71) 



Incidentally, we point out that the groundstate degeneracy always decreases under renormalization, which is the 
content of "g-theorem" [8]. This leads to a zero temperature entropy for the impurity given by Sj mp (0) = ^ In 2. 

Temperature corrections to this result have also been calculated [60] by using real fermion coordinates (Majorana 
fermions) ip a (x)(a = 1,2,3), which obey the anticommutation relations {ip a (x), ip b (y)} = S ab S(x — y), describing the 
relevant coordinates only. This corresponds as well to a regularization scheme where the fixed point has been moved 
to J = oo. The VT behavior is recovered, as metioncd in the introduction to this Section (see eq.(52)). 

C. Can we reach the two channel Kondo flxed point in a quantum dot ? 

A possible experimental realization of the two-channel Kondo fixed point in a quantum dot has been recently 
proposed [24]: exact diagonalization results of a vertical quantum dot with five electrons show that it can be tuned, 
by means of a strong external magnetic field, at the degeneracy point between energy levels with S = 1/2, but with 
different orbital angular momentum. Vertical cylindrical contacts provide single particle energy subbands labeled by 
the cross-sectional angular momentum and the k vector of the incoming/outgoing electron. Appropriate tuning of 
the electron density in the contacts offers the chance of including two electron channels only. The advantage of this 
setting is that no exchange coupling can take place between the channels due to symmetry, so that they act as totally 
independent. Selection rules due to the cylindrical symmetry enforce angular momentum m conservation and spin 
component a conservation in the cotunneling processes. For the special setting of the proposed device only a hole (h) 
process for J. —spin and a particle process for f —spin are allowed. They differ in the sign of the potential scattering 
but this difference is inessential. In fact, a particle- hole transformation on the fermion fields of the f channel only 
reverses the unwanted sign without affecting the exchange coupling. This puts both channels on an equal footing and 
points to an "orbital" Kondo coupling where the spin acts as the label for the channel. At this stage the dot plays 
the role of an effective spin 3/2 impurity, interacting with two channels. So, as it stands the system would flow to an 
undcrscreend situation in lowering the temperature. However, provided that Zccman spin splitting A is larger than 
the Kondo temperature for the underscreened fixed point Tj£, there is a crossover temperature T c at which two of 
the four S = 1/2 states are ruled out of the cotunneling processes, thus allowing for an effective spin 1/2 two channel 
Kondo flow. 

Detecting the effect requires an appropriate control of the hybridization of the dot with the contacts (so that 
A >> TJf ) and a proper tuning of the gate voltage V g . This allows for fine tuning of the exchange couplings between 
the dot and the two channels. It is well known that the hardest condition to be realized is total equivalence of the 
channels in the coupling. Were this not the case, the system would prefer one of the two channels as the dominant 
one and flow to a more conventional one-channel Kondo state. By tuning V g one can cross the point where the two 
channels are equally coupled which allows for scaling towards the two-channel fixed point. The measured conductance 
g (T) vs T will exhibit quadratic behavior at low temperature except for a crossover to a square root behavior, at 
the appropriate V g value which makes the two channels perfectly equivalent (sec Fig. 6). The delicate point in this 
proposal is that the requirement of full cylindrical symmetry of the device is essential. 

Another setup for measuring the two channel Kondo conductance has been proposed recently [25]. This proposal 
relies on an extra contact to provide the second channel. It is recognized however that this would introduce cross 
exchange terms between the three leads which do not allow for equivalent coupling of the two channels. In fact, any 
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diagonalization of the scattering problem of the kind of the one outlined in the previous Subsection to isolate the 
independent channels can never produce equal coupling as long as off-diagonal terms are non zero. To overcome this 
difficulty the authors propose that the extra contact is a larger dot itself, having a charging energy E c which hinders 
exchange coupling with the other two leads. Of course its size should be appropriately tuned. A smaller size implies 
a larger level spacing 6 inside it and S should be low enough because it acts as the low energy cutoff in the scaling 
toward the fixed point. Therefore, there is a delicate trade-off between incrcsing E c to prevent cross exchange terms 
and decreasing 5 not to stop the flow when reducing the temperature. 

We have stressed that the NFL fixed point in the two channel Kondo coupling at T — can only be reached by 
tuning the exchange couplings of both channels Ji, J2 exactly equal. Hopefully this can be done by changing an 
appropriate gate voltage across a critical point V* . According to Fig. 4 the systems flows to J\ — > 00, J 2 — > for 
V g < Vg* and to J2 — > 00, J\ — ► for V g > V g * , both of which are Fermi liquid fixed points. It has been argued 
that the quantum transition across V* should display a quantum critical region in the T — V g plane whose critical 
properties can be determined [26] . This offers better chances to spot whether we are in the vicinity of the NFL fixed 
point or not, also at finite temperature. 

VIII. SUMMARY 

In this review we focused on equilibrium transport properties across a quantum dot (QD). A QD is a tunneling 
center for electrons coming from the leads. Depending on the transparency of the barriers and on the temperature, 
the coupling V to the leads is weak or strong. 

In the weak coupling regime, tunneling can be dealt with perturbatively. In the case of the QD, a Coulomb Blockade 
zone is delimited by two sharp conduction peaks which grow as |V^| 2 /T in lowering the temperature. Conduction 
inbetween is due to cotunneling processes and is exponentially damped in temperature. Differential conductance is 
quite small, being ex |T^| 4 /J7, so that the charge degree of freedom is freezed on the dot. 

In lowering the temperature, non perturbative coupling of the dot to the delocalized electrons of the contacts can 
occur and conductance can reach the unitary limit in a CB valley. For pedagogical reasons we have reviewed the 
old-fashoncd Anderson- Yuval model for the correlated state. Otherwise we have used the poor man's scaling approach 
in the regime where perturbative scaling applies and the Nozicres scattering approach which entails the fixed point 
physics at zero temperature. These approaches do not allow for quantitative results which are better obtained with 
the numerical rcnormalization group (NRG), real time RG and Non Crossing Approximation (NCA) (the last one 
preferably in the overscreened case when a non Fermi liquid (NFL) fixed point is reached), but they offer a more 
transparent view of what is going on. 

We have briefly reviewed the various types of Kondo couplings. 

We have considered the case in which the GS of the dot is degenerate because of spin: if temperature is low enough 
(T < Tk, with Tk depending on the transparency of the barriers), spin flip processes proliferate and the magnetic 
moment of the dot is partially or fully screened . In this case an applied magnetic field has a disruptive effect on the 
Kondo peak of the conductance. The interesting case of a crossing between different dot spin states induced by the 
magnetic field has been also discussed. 

We have also reported on other possible realizations of Kondo physics involving orbital degrees of freedom (orbital 
Kondo). The most favourable set up for this case is a vertical geometry of the dot and the contacts with cylindrical 
symmetry. In this geometry a magnetic field orthogonal to the dot (which may be strong) can induce level crossings 
and produce the degeneracies of the dot GS which is required for Kondo conductance to take place. Hence one can 
have Kondo effect also with an even number of electrons on the dot and zero total spin. 

Some attention has been devoted to the multichannel Kondo effect. The overscreening case can lead to a NFL 
fixed point at zero temperature. In particular we have discussed the two channel spin 1/2 Kondo state and reported 
on possible experimental realizations that have been proposed. It emerges that conditions to be met are very tough. 
Nonetheless its achievement would probe a beautiful piece of the physics of the strongly correlated systems. 

IX. APPENDIX A: RG EQUATIONS FOR THE COULOMB GAS MODEL OF EQ. (29) 

We summarize here the scaling of the action in eq. (29) in order to find out the behaviour of the system at large time 
scales (low temperature) [21] [34]. The cutoff is rescaled according to: tm — * (1 + A)tm with A = Atm/tm ~^ din tm- 
This adds a factor e - 2N H 1 - a / 2 ) xhc first term arises from rf^ in the denominator, while the one oc a 2 arises from 
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the In— term in the interaction term. Indeed charge neutrality implies that = (X^<7*) 2 = ^2ijtj QiQj + Si Qi an d 
qij = ±1. This factor rcnormalizes the fugacity Y — > y + dy, with: 



y + dy «y e A(1 - a2/2) , 



(72) 



what gives the first of eqs. (30). The interaction strength a is renormalizcd by flip-antiflip (particle-antiparticle) 
fusion. 

Let us consider now all the configurations in which pairs of neighboring charges qi at Tj and qj — — qi at Tj (where 
j = i ± 1) are at a distance between tm and tm(1 + A). In increasing the scale, these pairs are seen as a neutral 
compound which screens the interaction between other charges ("fusion " of pairs). Let us consider one single fusion 
process and develop that part of the action that contains their coordinates: 



. 5 (2iV + 2) = _ 5 (2iV) + ^. lll 
hi 






tm 
Tk-r- e/2 



tm 



-In 



T k -T + e/2 



tm 



(73) 



Here we have defined e/2 = T ' 2 Tj and r = T> 2 TJ ■ The integral over r and e for e/rj| << 1 which appears in the 
partition function of eq. (29) is: 



? (2JV) 2 /^h- 1 -™ dr /-™ (1+A) de 
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where the term with k = i — 1 (fc = j + 1) in the first (second) sum vanishes. Doing the same for each interval i — j 
and summing, each term is repeated twice. This generates an extra contribution to the action that can be interpreted 
as the renormalization of the coupling constant: 



a 2 -► a 2 + da 2 ; da 2 = -2XY 2 a 2 . 
Because A « d In tm, the second one of eqs. (30) is obtained. 



(75) 



X. APPENDIX B: TEMPERATURE DEPENDENCE OF THE CONDUCTANCE 

Ohm's law can be derived semiclassically from a Boltzmann equation which accounts for weak scattering of the 
charge carriers against impurities and defects, when driven by an electric field. Scattering provides a mechanism for 
relaxation from a non equilibrium to a steady state flow [62] . We follow here a simplified approach starting from the 
velocity of band electrons v k — Ve^/ft.. The current density is: 



J 



2e 
3V 



X>(/*-/ fc ) 



where fk(fk) is the non-equilibrium (equilibrium) Fermi distribution and V is the volume. Now we take: 



fk - fk 



dj°_ 

de k 



Sek 5e k = -eE-v k T(k), 



which defines the relaxation time r(fc), so that, according to Ohm's law, we get: 

° = -wY k - VkT{k) de- k - 



(76) 



(77) 



(78) 



In th case of s— wave scattering a spherical wave is diffracted from the impurity with maximum amplitude; this 
increases the flux propagating backward. Assuming that quantities depend on e only and are mostly evaluated at the 
Fermi level we have (T <~ 0) : 
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2e^ 
3V 



v(0)v% f der(e,T) 



df°_ 
de ' 



(79) 



In this way we recover the simple microscopic formula for Ohm's conductivity, first derived by Drudc in 1900, 
a = ne 2 r /to, where r is an average relaxation time. The latter can be defined in terms of the mean free path 
between two scattering events I = vft or distance between impurities. This formula is valid in the limit k F l >> 1 
which we assume to be the case. In such a limit no localization effects occur and one can show that it is dp/dT > 
always, which is violated in the case of Kondo conductivity in diluted alloys. 

In fact, for the large U Anderson model, the explicit dependence of the transport time on temperature at the 
Fermi energy has to be taken into account separately. In addition, the Sommcrfeld expansion can be used at finite 
temperatures [10]: 



dc 



d.f{e) 
de 
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(80) 



The main ingredients to derive this formula are the Fermi liquid nature of the zero temperature GS and the Taylor 
expansion of the temperature dependence of the relaxation time. The prefactors have been specialized to the case of 
the symmetric Anderson model (so that the lifetime of the Kondo resonance is T = 4fcsT/f/7r). 
It follows that, in diluted alloys, the resistivity takes a maximum at T = 0: 



P(T) 
P(0) 



1 



ir 4 T 2 
167f 



(81) 



The extra correction on the r.h.s. of eq. (80) arises from inelastic scattering at finite temperature, as shown by 
Nozieres [40], according to the following argument. 



In an elastic resonant scattering process at T = is — ^sm{t°} 
finite T, we have to consider the total relaxation rate, defined as: 



v(0)tT 



sin 5°. In a diffusive s— wave scattering at 
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rje 



—TT = 7 ,Wfcfc'(l 



k 1 



2ir 

k-k') = — NM0) 

n 



de k ,5{e k -e k ,)\ < k\t\k' > | : 



(s— wave scattering implies that the correction k-k' averages to zero). The inelasticity can be accounted for by defining 
an effective elastic S— matrix which is now no longer unitary and an effective phase shift which is different from the 
one at T = 0. 



The full S— matrix, instead, is of course unitary: 1 
27riv(O)t* a0 8(e a ~ e p )] = E/3^/3 



E^ 2 {v (Q)f \t aft \*8{e a 



JlelSapl 2 = J2 [8a/3 ~ 2Trw(0)t a gS(e a - e )][S a p + 



ep). Now, we separate the elastic channels (3 = a 



from the inelastic ones on the r.h.s. and we introduce the inelastic cross-section W™ = 2^^^ , Q f(0)\t a p\ 5(e a — £/3); 
hence, from unitarity it follows that 1 + 4tt 2 (^(0)) \t aa \ 2 + 2irv(0)W^ 1 = 1. Thus, we can define an effective phase 
shift (which we call S again) in presence of inelastic scattering by introducing an effective "elastic" S matrix: 



in\ 2 J2i8 



1 - 2mu(0)t aa = (1 - 2tw(0)W") 

where the square root can be expanded. 

According to the optical theorem, we can write the total transmission probability as: 



W a 



-2%m{t aa } = W™ cos 25 + 



(1- cos 26) 

i/(0)7T 



(82) 



(83) 



Now, let us notice that, because S ~ 7r/2, we have cos 25 < 0. Therefore, for each channel a the quantity W a is always 
smaller than the elastic case value: 2 sin 6/(v(0)ir). But one can thermally average over the a channels; according 
to the usual phase space argument which is invoked when scattering occurs close to the Fermi surface at very low T, 
we obtain an average W m proportional to T 2 : W m ~ AT 2 . This gives rise to the first term in the expansion of eq. 
(80). The second term can be viewed as contributing to the Sommcrfeld expansion with 5(e) = ^tt — ae from which 
it follows that sin 2 6 ~ 1 — a 2 e 2 . 

In the case of resonant tunneling across the dot, the conductance is given by eq. (20). Using the result of eq. (83) 
directly one obtains that conductance has a maximum at zero temperature and first corrections in temperature are 
again of 0(T 2 ) (see eq. (49)). 
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Figure Captions 

• Figure 1: Schematic drawing for a dot and the contacts in a vertical setup; possibly a magnetic field is applied 
along the axis. 

• Figure 2: Differential conductance on a gray scale as a function of both V g and Vds] the Kondo effect shows up 
near V ds = [29]. 

• Figure 3: Renormalization-group flow diagram for small J [21]. 

• Figure 4: Qualitative renormalization-group flow diagram for the anisotropic two channel Kondo problem; the 
non trivial fixed point corresponds to J\ = J2 = J* (channel symmetry) [42]. 

• Figure 5: Plot of the universal function / (x) vsi = T /Tk [11]. 



• 



Figure 6: Sketch of the conductance across the dot as a function of the temperature T and the gate voltage V g 

[24]. 



• Figure 7: Second-order diagrams T>\ and T>2- Conduction fermions are represented as full lines, while dashed 
lines represent propagation of d- fermions (dot's states). 

• Figure 8: Third-order vertex renormalization: the first one (P) is a "parquet-type" diagram. Its contribution 
is accounted for in the integration of the second-order RG equation. The second one (NP) is a "non parquet" 
diagram. It provides an additional third-order contribution to the /3-function. 
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